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I. INTRODUCTION 

Understanding strong quantum fluctuations continues 
to be a formidable challenge. Where two (or more) 
phases compete and are separated by a continuous quan- 
tum phase transitioi>i, i.e., at zero temperature, such 
fluctuations are particularly strong. Generically, such a 
quantum phase transition is signaled by the decay of ele- 
mentary excitations. Far away from the phase transition, 
spectroscopic probes show dominant sharp (5-peaks at low 
energies which result from stable elementary excitations, 
so-called quasi-particlcs. But on approaching the phase 
transition, the spectral weight in the dominant quasi- 
particle peak is reduced further and further and shifted 
to contributions of more quasi-particle. Multi-particle 
spectra with considerable weight are an important signa- 
ture of dominant quantum fluctuations in general. The 
vanishing of the single quasi-particle peaks at zero tem- 
perature is a smoking gun for a quantum phase transition 
in particular. 

In this context, the transverse field Ising model 
(TFIM)^ is a popular generic model describing mag- 
netic excitations and displaying a quantum phase tran- 
sition between the disordered phase in the limit of dom- 
inating transverse field and the ordered phase in the 
limit of dominating longitudinal Ising coupling. Due 
to its relative simplicity, the TFIM provides a con- 
venient test case in the development of theoretical 
approaches^. This is particularly true for the one di- 
mensional case, for which fcrmionization by a Jordan- 
Wigner transformation^ yields an exact solution^"—. 

The calculation of dynamical correlations in the TFIM 
is an active field of research. While transverse correla- 
tions can be treated in terms of fermions^i^, longitudinal 
correlations require different approaches, because of their 
non-locality in the fermionic picture. Based on an equa- 
tion of motion for the longitudinal correlations^^ there 
has been a series of papers inve stig ating the scaling region 
around the critical point Ref. [lll - llO . In 2009, Perk and 



Au-Yang computed results for the time-dependent longi- 
tudinal correlation functions by solving the coupled dif- 
ferential equations and complementing them with long- 
time asymptoticsiii. But so far no momentum and fre- 
quency resolved analysis has been performed, which also 
applies away from the scaling regime. 

Quantum magnets in the vicinity of quantum phase 
transitions are dominated by strong quantum fluctua- 
tions. Quantum fluctuations are strongly favored by 
low-dimensionality and by frustration. Theoretically, a 
particularly clear sign of dominant quantum fluctuations 
is the fractionalization of elementary excitations. For 
instance, conventional spin waves (magnons) split into 
two spinous in antifcrromagnetic Heisenberg chainsi^Ti^i. 
Before complete fractionalization occurs, the spectral 
weight, as observed in inelastic neutron scattering, shifts 
from the channel of a single elementary excitation to the 
channel where two and more elementary excitations are 
crcated^^i^. Thus, also the experimental focus is di- 
rected more and more to continua formed bymore than 
one excitation, see for instance Refs. [TqI. [20l l24l and[25l. 

In view of the above considerations, the present arti- 
cle pursues two goals in a study of the one-dimensional 
(ID) TFIM. First, we show how the special algebraic 
structure ('string algebra') of the operators occurring in 
the ID TFIM enables its solution by a continuous uni- 
tary transformation (CUT) in the high-field phase to very 
high accuracy. Upon completion of our calculations, we 
learned that this algebra was introduced and used before 
m an algebraic solution of the TFIM^^. 

This algebra paves the way to treat a larger class of 
models of which the Hamilton operators can be expressed 
by operators belonging to the string algebra, for instance 
XY models in transverse fields. Second, we compute 
the three-particle contributions to the diagonal dynamic 
structure factors (DSF) in the CUT framework in the 
high-field phase. To our knowledge, this is the first time 
that these subdominant contributions are computed, ex- 
cept in the scaling region around the quantum phase 
transition. Thereby, interesting predictions for future ex- 
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perimental studies are provided. 

In Sect. II, the model and known exact results are re- 
called. In the following section, the continuous unitary 
transformations (CUTs) are briefly reviewed. Sect. IV is 
devoted to the algebra of string operators which are sub- 
sequently employed. Sects. V and VI comprise our static 
and dynamic results, respectively, while the conclusions 
are drawn in Sect. VII. 



II. MODEL AND EXACT RESULTS 

The Hamiltonian for the transverse field Ising model 
(TFIM) reads 



TFIM 



where the cr" arc the Pauli matrices and the sum i runs 
over all lattice sites. We normalized the distance between 
two sites to one. The transverse field strength is given 
by the parameter F while J denotes the strength of the 
antiferromagnetic coupling between two adjacent sites. 
The antiferromagnetic exchange can be converted to a 
ferromagnetic exchange J — — J by a tt rotation around 
for every second site i. This translates to a shift of tt 
in momentum space. 

The model has a quantum phase transition at J = 2r 
and it is self duali. Similar to Ref. [l], we introduce the 
parameter x = J /2T. The starting point for the CUT 
calculations is J = 0. Hence, an elementary excitation is 
given by a single spin flip. For finite J the energy of these 
excitations becomes momentum dependent. We will refer 
to these excitations as quasi-particles. We expect that 
the perturbative ansatz breaks down once we reach the 
critical value J = 2F. Hence we focus on the static and 
dynamic properties for J < 2F. 

The model was solved exactly by Pfeuty in 
1970i, based on the work by Licb ct ali^ and 
Niemeijer— . Pfeuty's solution uses the Jordan- Wigner 
transformation^ to map the Hamiltonian in Eq. ([T]) to a 
chain of free fermions which is diagonalized by a Bogoli- 
ubov transformation^^. This approach yields the exact 
expression for the ground state energy per site 



|l = / ..(,)dg. 



where denotes the dimensionless dispersion 



j{q) = + x'^ — 2x cos(g) . 



(2) 



(3) 



The dispersion with dimension is given by Tuj{q). 

From the dispersion we can easily extract the energy 
gap of the lowest lying excitations 



(4) 



Another interesting quantity worked out by Pfeuty is the 
transverse magnetization 
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where \g) denotes the ground state of the TFIM. In the 
following sections we will compare our results with these 
exact expressions in order to validate the CUT approach. 

Beside the static properties stated above, dynamic 
properties are important in order to explain experimental 
results. Although the TFIM is analytically integrable, 
the evaluation of longitudinal dynamic correlations re- 
mains a very difficult task which requires considerable 
numerics, see for instance Refs.ll7lor[29l In the fermionic 
picture, this is due to the non-locality of the Jordan- 
Wigner transformation. 

One important quantity in the study of spin systems 
is the dynamic structure factor (DSF) 



S"''{uj,Q) 
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where a, /3 S {x, y, z}. Here Q denotes the total momen- 
tum and uj the frequency. The DSF is directly linked to 
the differential cross section in inelastic scattering exper- 
iments, see for instance Ref. [s^. 

Due to the symmetry af -> —erf of i/TFiM no correla- 
tions occur for a = x, P = z and a = y, (3 = z and vice 
versa, but for a ^ x, (3 ~ y and for a = /3 the DSF may 
and will obtain finite values. 

In the following, wc focus on the diagonal part of the 
DSF, i.e., a = p. For a = z exact expressions are 
know n^i^i^^ , because the observable cr^ remains local in 
the Jordan- Wigner representation of the TFIM. At zero 
temperature case, this expression reads 



S'^Q,uj) = J dfci [l-/(Q,fci)] 

— TT 

<5(w~w(fci-Q/2)-w(fci+g/2)), (7a) 



with 



(F + ^ cos(fci - Q/2)) (F + f cos(fci + Q/2)) 



uj{ki - Q/2)iu{ki + Q/2) 



(7b) 



It consists of a spectral density of scattering states of two 
elementary excitations with total momentum Q. For a = 
X and a = y, only the one-particle contributions have 
been calculated by Hamer et al. in 2006^. They used 
series expansion techniques to propose the expressions 



sriQ) 
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5r(g) = [i-x^l^^Q) 
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for the one-particle contribution to the equal-time struc- 
ture factor. By comparing their results to correlation 
functions in the two-dimensional classical Ising model, 
see Ref . [H and [13, they could show that the expressions 
above are indeed exact. Hence the full one-particle struc- 
ture factor is given by 



(9) 



For higher-particle contributions to the longitudinal 
DSF much less is known. In 1978, Vaidya and Tracy 
computed exact expressions for the longitudinal correla- 
tion functions in the anisotropic XY model in the time 
domain^^. They evaluated the resulting expression in 
frequency space up to the three-particle contributions. 
But their results are limited to the scaling region at low 
energies, very close to the critical point. Furthermore 
Miiller and Shrock calculated frequency-integrated wave 
number dependent susceptibilities for the TFIM at the 
critical point in Refs. [Tsl and [l^. Our results will be 
complementary to theirs. 



III. CONTINUOUS UNITARY 
TRANSFORMATIONS 

We use the method of continuous unitary transforma- 
tions (CUT) to derive effective models which allow for an 
easier evaluation of static ground state properties and dy- 
namic correlation functions. The idea of CUT was first 
introduced by Wegner— and independently by Glazck 
and Wilson^i^a. 

The concept of CUT is to systematically finda uni- 
tary transformation that maps the Hamiltonian to a di- 
agonal representation. One introduces a family of uni- 
tary transformations U{1) depending differentiably on 
a parameter I G IR"^. The unitary transformation is 
characterized by its anti-Hermitian generator ri{l) ~ 
{diU Then, a short calculation yields 
the flow equation 



diH{l) = [r,{l),H{l)] 



(10) 



for the /-dependent Hamiltonian H{1). In general, it rep- 
resents a system of coupled differential equations for the 
prefactors of all operators appearing in H{1). We refer 
to it as the differential equation system (DES). Without 
further truncation, the DES generically comprises an infi- 
nite number of variables. In practice, various truncation 
schemes help to keep the DES finite. For / — )■ oo the 
Hamiltonian acquires its final form and it is denoted as 
the effective Hamiltonian 



H, 



off' 



H{1)\ 



C/(oo)iJC/t(oo). 



(11) 



The convergence for £ — >■ oo is assumed; it cannot be 
proven generally for infinite dimensional quantum sys- 
tems because it depends on the specific form of the gen- 
erator as well as on the employed truncation scheme. 



Note that observables O also need to be transformed to 
effective observables by the same unitary transformation. 
This results in the flow equation for observables 



diO{i)^m,o{i)] 
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which yields the effective observable Oofi for Z — > oo. 

The generator characterizes the CUT and the flow of 
the Hamiltonian. Thus, the choice of the generator is 
an important issue and it still represents an active field 
of research, cf. Refs. 0, IH, l39l - l4a In this paper we use 
the (quasi-) particle conserving (pc) generator, which was 
first proposed by Mielke^ in the context of banded matri- 
ces and independently by Knetter and Uhrig^ for many- 
body problems. By 'quasi-particle' we mean the elemen- 
tary excitation. The pc generator directly aims at these 
quasi-particles. The goal is to eliminate terms that do 
not conserve the number Q of quasi-particles 



0. 



(13) 



The pc generator is given in matrix representation in the 
eigenbasis of Q by 



VpcAjil) = sgn{qi - qj)hij{l), 



(14) 



where qi denotes the eigenvalues of the operator Q. 

An equivalent description of the pc generator can be 
given by decomposing the Hamiltonian into parts that 
create, H^{1), conserve, H^{1), and annihilate, H^{1), 
quasi-particles. Then the Hamilonian reads 

H{l) = H+{l) + H"{l)+H-{l) (15) 

and the quasi-particle conserving generator 



7]pc=H+{l)-H-{l). 



(16) 



The convergence of the flow induced by this generator is 
proven for finite-dimensional systems; extensions to infi- 
nite systems are also available^. The pc generator pre- 
serves the blockband diagonal structure, i.e., the maxi- 
mum number of particles created or annihilated docs not 
change during the fiow^Si^i^. 

The CUT method consists of two basic steps. The com- 
mutator in Eq. ([TU]) needs to be calculated, followed by 
the integration of the resulting flow equation. The latter 
can easily be done with standard numerical integration 
algorithms or even analytically. 

In general, commuting H with rj creates new types 
of terms which were originally not part of the Hamilto- 
nian. For systems in the thermodynamic limit, all sorts 
of new terms may arise connecting more and more sites 
over larger and larger distances. In a numerical calcu- 
lation we cannot treat an infinite number of operators, 
hence we have to restrict ourselves to operators which 
are physically relevant. In this paper we use the previ- 
ously introduced directly evaluated enhanced perturbative 
CUT (deepCUT)i^. The idea of deepCUT is to trun- 
cate operators and contributions to the DES according 
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to their effects in powers of a small expansion parameter 
X. Roughly speaking, the order n in x is the truncation 
criterion. More precisely, a certain contribution to the 
DES is kept if it affects the targeted quantities (here: 
ground state energy and one-particle dispersion) in order 
m < n'va X. Details can be found in Ref. Ilsl 

Thus we write our initial Hamiltonian in the form 

H = Ho + xV (17) 

where Hq describes the unperturbed Hamiltonian and V 
represents a perturbation. We expand the operators in 
the basis {Ai} which is chosen such that the effective 
Hamiltonian can be computed exactly^ up to order n in 
the parameter x. Then the flowing Hamiltonian can be 
denoted as 

H{l) = Y,h^{l)A, (18) 

i 

where the prefactors hi{l) depend on the flow parameter 
I. For the generator we choose the same operator basis 
with the same prefactors 

i i 

where rj [■] is a superoperator applying the generator 
scheme. For the pc generator, rj [Ai] = Ai if Ai creates 
more quasiparticles than it annihilates, [Ai] = —Ai if 
Ai annihilates more quasiparticles than it creates, other- 
wise [^i] = 0, cf. Eq. p^ . 

With this definitions, we obtain the flow equations 

dMl)^Y.^^Jkhj{l)hk{l). (20) 
j,k 

We call the Dijk E C the contributions to the DES. They 
are obtained in a perturbative calculation up to order n 
by calculating the commutator in Eq. (jlOp and expanding 
the results in the chosen operator basis. Note that the 
numerically evaluated DES also comprises powers in x 
beyond the order n^. 



IV. STRING OPERATORS 

In the previous section, we explained how continuous 
unitary transformations are applied in a general context. 
Here we specify the approach for the transverse field Ising 
model. For the TFIM in the high-field limit, we use the 
the state with all spins down 

|o) = i---;j-i;,;,+i---), (21) 

as the reference state, i.e., as the vacuum of elementary 
excitations. This corresponds to the strong field limit 
r 00 in the TFIM, which is also the starting point for a 
perturbative approach in the parameter x = J/{2T). An 



elementary excitations, i.e., a quasi-particle, is created 
by the spin- flip operator CT;+. We denote this state by 

|0=a+|0). (22) 

It is obvious that no two excitations can be present at the 
same site so that the quasi-particles behave like hardcore 
bosons. But multi-particle states can straightforwardly 
be created by flipping spins at different sites. 

These ideas suggest a basis of operators consisting of 
monomials made from the local operators 

K'-^.-^-^.-l}' (23) 

where stands for particle creation, a~ for particle an- 
nihilation, and o'j'crj counts whether a particle is present 
at site j {(T^crJ = 1) or not (crfaj — 0). This approach 
is in line with the general structure explained in Ref. 1461 : 
we call it henceforth the multi-particle representation. 

The number of such monomials grows exponentially 
with the number of sites which are non-trivially involved 
because at any site a quasi-particle may be created or 
annihilated or simply counted. For each additional site 
occurring in the course of the iterated commutations, the 
number of operators to be tracked grows roughly by a fac- 
tor of 4 (neglecting reductions due to symmetry effects) . 
This is a major drawback if one aims at higher orders. 
Therefore, we introduce a simpler modifled operator ba- 
sis, which we call string algebra^ which is more appropri- 
ate for the TFIM, see also Ref. [H. We stress that the 
possibility to introduce a string algebra is connected to 
the Jordan- Wigner representation of the Hamiltonian in 
terms of non-interacting fermions. 

The string algebra consists of operators which are given 
by the following product of Pauli operators 

(j+n-l \ 
n < (24a) 

3 

with {0, e} G { + !^} and n € N. Each string operator 
consists of a product of adjacent operators, framed 
by spin flip creation- and/or annihilation operators. The 
Gz operators form the string between the pair of spin 
flip operators. We refer to n as the spatial range of an 
operator. In contrast to the local set of operators used in 
the multi-particle representation (|23p the string algebra 
is more transparently expressed by the set 

{a+,a7,a|,l}. (25) 

The key point of the string algebra is that excitations 
or annihilations of quasi-particles occur only at the end 
points of the string. Thus, for given end points, there 
are only four string operators to be considered. If exci- 
tations or annihilation could occur anywhere along the 
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string one would have exponential growth of the number 
of operators. 

In Eq. (|24|) . we defined the translationally invariant 
form of string operators. When dealing with loeal observ- 
ables, it is also useful to introduce local string operators 



O 



' j+n-1 

.k=j+i 



(26a) 



(26b) 



with {(/>, e} G {+7^} and n G N . Note that a transla- 
tionally invariant string operator is given by the sum of 
local string operators. 

It is also useful to define a string operator of range 
consisting of a single cr^ matrix. 



Tn 



(27a) 
(27b) 



For n — 1, the definitions and (|26p correspond to 
a normal hopping term or pair creation or annihilation 
operator. These operators cannot be distinguished from 
operators in the multi-particle representation^. 

For n > 1, the situation is different. For example, in 
the case n = 2, (jj = + and e = — can be re-expressed in 
multi-particle representation by 

T2^- = E<'^l+i^7+2 (28a) 
j 

j 

= E (2«+i'^7+i'^7+2 - 4^7+2) ' (28c) 
j 

which is the sum of a quartic interaction term and a hop- 
ping term, because we re-expressed a|^]^ = 2cr^j^cr^-^ — 1. 

This simple example illustrates the computational ad- 
vantage of the string algebra. If we tracked all opera- 
tors in multi-particle representation, a single string oper- 
ator of range n would require 2"~^ multi-particle opera- 
tors, clarifying the previous statement on the exponential 
growth of the number of such monomials. Therefore, if 
a model can be diagonalized within the string algebra, it 
is highly advantageous to describe all operators in terms 
of string operators. 

With the above definitions, the Hamiltonian of the 
transverse field Ising model is formulated in terms of 
string operators 



i^TFIM = 2 E"^! + fE + '^t'^t+l + ^-C- 



J, 

A' 

J J 



^To + ^ (T+- + T-+ + T++ + T- 



(29a) 
(29b) 



Next, we study the action of hopping terms on single 
particle-states 



T-+ \l) = Y.''J^Ui^h2 ■ ■ ■ 4+n 10 (30a) 

3 

= E^'-^-^I+i^I+2 • • ■ 4+n |0) (30b) 

3 

= {-1)"-^ \l + n) . (30c) 

In the second line we used the property \l) = (5;.j |0) 



and we know that <t| |0) 



which yields the final 



result. If there is only one quasi-particle in the system, 
the only difference to conventional hopping is the factor 



(-1)' 



For subspaces with more quasi-particles, we 



have to take into account that there may be particles 
on the sites l + l,l + 2...l + n— 1. They modify the 
exponent of (—1) and can thus change the sign of the 
resulting state. 

In order to apply the deepCUT to the TFIM, we have 
to calculate the contributions to the DES in Eq. ([20]) . 
Thus, we calculate the commutator between two opera- 
tors of the Hamiltonian H and the generator 77. In Ap- 
pendix [X] we show that the string algebra is closed under 
such commutations. This means that the commutator 
of two string operators can again be written as a lin- 
ear combination of string operators. The closure of the 
string algebra allows us to set up the flow equation in 
very high order because the number of operators to be 
tracked grows only linearly for the Hamiltonian. For local 
observables within the string algebra it grows quadrati- 
cally which is still a moderate growth. This is the key 
observation of the present article. 

Explicitly calculating all distinct commutators of 
string operators, see Appendix |X1 allows us to determine 
all contributions to the DES analytically. In this way, we 
calculate the flow equation up to infinite order in x. We 
parametrize the Hamiltonian 



i?TFIM(0 = to{l)n + E 4^(0 {T+- + h.C.) 
n=l 

00 

+ Y.t++{l){T:++h.c.), (31) 

n=l 

and the generator for the CUT 



(32) 



In Appendix IB] we derive the flow equation for the pref- 
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actors io,i^ i^^"*"- It reads 

oo 

dit, =2^(4+)', (33a) 

n=l 

A;+i=m \k-l\=m 

da+-^2 J2 *r^r-2 ^ (33b) 

|/c— /| 

k-\-l—ra 

with m G N. Note that this is a differential equation with 
an infinite number of variables which grows, however, 
only linearly in the spatial range. This result is remark- 
able considering the fact that it would require tremen- 
dously more flow parameters if we formulated the prob- 
lem in multi-particle representation. Thus the string al- 
gebra allows us to evaluate the Hamiltonian transforma- 
tion up to very high orders, which is especially important 
on approaching the quantum critical point x = J/2T = 1. 

V. STATIC RESULTS 

In this section we evaluate and present static results 
for the transverse field Ising model. The expression 
'static' refers to time-independent properties. We treat 
the ground state energy per site in Sect. IV Al the mag- 
netization in Sect. IV Bl and the momentum-integrated 
spectral weight in Sect. IV Cl aud the momentum-resolved 
static structure factor in Sect. IV Dl 

A. Ground state energy 

Due to the only linearly growing number of string op- 
erators, we are able to obtain the ground state energy 
per site up to order 256. Higher orders do not improve 
the results significantly so that we restrict ourselves to 
orders up to 256. 

Figure [T] compares the exact result for the ground state 
energy per site to CUT results in various orders in x. 
As expected, the accuracy increases for increasing order. 
Even close to the critical point the CUT result of order 
128 and the exact results can barely be separated. The 
inset shows the deviations from the exact results. On 
the logarithmic scale, straight lines indicate power laws 
for these deviations as expected in a perturbatively con- 
trolled approach. We checked that the slopes of the lines 
correspond to the order of calculation by fitting the de- 
viations to power laws, see Tab. ID 

From the inset we also read off that the decpCUT in 
order 128 calculates the ground state energy per site cor- 
rectly to the fifth digit, even at the critical point. The 



Order 


Exponent 


Fitting Error 


9 


10 


± 3 


16 


21 


± 4 


32 


33 


± 3 


64 


72 


± 6 


128 


132 


± 5 



Table I. Exponents of the power laws for the deviation of the 
ground state energy obtained by numerical fits. 
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Figure 1. (Color online) Ground state energy per site as func- 
tion of J. Comparison of tire exact result to CUT results in 
various orders of a; = J/(2r). The inset shows the abso- 
lut diflference between the exact result and the CUT calcula- 
tion on a logarithmic scale. The critical point is located at 
logio iw) = 0. 



calculation in order 256 is not shown in the graphs be- 
cause it would be indistinguishable from the other curves. 
It improves the result in order 128 at the critical point 
by about one digit. 



B. Transverse magnetization 

Next, we examine the transverse magnetization 

M^--^J2(9\^^\9)- (34) 

j 

Here \g) denotes the ground state of the Hamiltonian 
which is mapped to the zero-particle state by the CUT. 
Note that the expression 'transverse' refers to the direc- 
tion of the external field, which is the z-axis in our model. 
In the limit J — >■ 0, all spins are aligned along the exter- 
nal field, so that = 1 holds. 

For J > 0, the spins are perturbed by the antiferro- 
magnetic interaction which reduces the transverse mag- 
netization. To obtain M^, we transform the observable 
cr| by the continuous unitary transformation to an effec- 
tive observable. The operator of the transverse magneti- 
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Figure 2. (Color Online) Transverse magnetization as a func- 
tion of J. Comparison of the exact result with the CUT 
calculation. The inset shows a close view on the critical point 

J = 2r. 

zation can be expressed by a string operator 

^a|a = 0)=To. (35) 
j 

Due to this identity and the fact that the string algebra 
is closed under commutation, we know that the final ef- 
fective observable can be written as a linear combination 
of string operators 

a^(0 = 00(0^0 + + h.c.) (36a) 

n 

+ o++(0 (r++ + h.c.)]. 

None of the coefficients and contribute to the 
vacuum expectation value. But they can not be omit- 
ted during the flow of the observable. The transverse 
magnetization after the CUT reads 

M, = oo(oo). (37) 

Due to this simple form of the observable very high or- 
ders can be reached again. The transverse magnetization 
calculated by the CUT is compared to the exact result in 
Fig. [5J As expected the results improve upon increasing 
order. The largest error occurs at the critical point where 
the transverse magnetization displays a singularity. 

C. Spectral weight 

In this section, we discuss the CUT results for the 
momentum-integrated quasi-particle weight in the two 
diagonal channels and 5^^. We use the CUT frame- 
work to calculate effective observables which renders an 
easy evaluation for the spectral weight possible in vari- 
ous quasi-particle channels, see also Ref. [13. The total 



spectral weight can be split according to 

^"" = ^E(^rO=A+/2 + /3 + ..., (38) 

where /„ denotes the weight in the channel with n quasi- 
particles in the system. Introducing the CUT framework 
results in 

/r = (oi<eff<iffio>, (39) 

where cr" denotes the part of the effective observable 
which annihilates n quasi-particles. Since and cr^ cre- 
ate an odd number of spin flips, i.e., quasi-particles, and 
since the generator of the CUT preserves the parity of 
an observable, S*^^ and S^^ consist of 1,3,5,... quasi- 
particle contributions. With the help of the sum rule 
jaa = \^ can also check if our results are still valid 
for large values of x. 

We emphasize that the local observables af and 
transform into non-local operators under the Jordan- 
Wigner transformation which act on an infinite number 
of sites. Therefore, no easy analysis of these observables 
is possible in fermionic terms. An explicit evaluation re- 
quires either analytical mappings which enable an eval- 
uation in the scaling region^ or extensive numerics in 
terms of Pfaffians whose dimensions grow linearly with 
the spatial range of the correlation^^. 

The problem of an inflnite number of operators is 
avoided in the string operator basis ([25]) . But the calcula- 
tion remains cumbersome because the observables are not 
part of the string algebra and thus the structure of the 
effective observables is more complicated. This compli- 
cated structure prevents us from achieving very high or- 
ders, because the number of representatives to be tracked 
grows exponentially on increasing order. We are able to 
obtain results up to order 38. Then the computational 
effort reaches its limit in the present implementation be- 
cause the contributions to the differential equation sys- 
tem take more than 8 GB of memory and the number of 
operators to track is larger than 7 million. 

First, we address S^^ as function of J for which results 
are depicted in Fig. |31 The CUT results arc compared 
to the exact results from Ref. |H. The one-particle con- 
tribution shows a very sharp drop for J — > 2r with a 
singularity at the critical point. The CUT agrees very 
well with the exact results as long as the order of cal- 
culation is below the correlation length. Recall that the 
order is proportional to the range of the physical pro- 
cesses included in the calculation. For the calculation of 
effective observables within the string algebra this was no 
problem because large orders > 100 could be achieved. 
But for the longitudinal correlations, we obtain only or- 
der 38 so that particularly sharp edges such as the one 
in the one-particle spectral weight are not captured. But 
the agreement improves on increasing order. 

The spectral weight of the three-particle channel in- 
creases on approaching the critical point. Hence, spec- 
tral weight is transferred from the one-particle channel 
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Figure 3. (Color Online) One-particle and three-particle spec- 
tral weight as function of the parameter J. Comparison of 
the exact expression for the one-particle weight H8a|) with the 
CUT results. 



Figure 4. (Color online) One-particle and three-particle spec- 
tral weight as function of the parameter J. Comparison of 
the exact expression for the one-particle weight (|8b|) with the 
CUT results. 



to higher quasi-particle channels. The one- and three- 
particle terms are the dominant contributions to the total 
spectral weight for the parameters investigated. But for 
J > 1.9, the sum rule starts being violated in the CUT 
calculation. We attribute this violation to the calcula- 
tion in finite order. It appears that the CUT calculation 
overestimates the one-particle contributions close to the 
critical point. 

Next, we investigate S^"^ . This correlation is depicted 
in Fig. m in comparison to the exact result. Again, the 
one-particle contributions also vanish for J — >■ 2r. But 
the edge at the critical point is by far not as sharp as in 
S^^ because more spectral weight is transferred to higher 
quasi-particle spaces for lower parameters J. The sum 
rule is again violated for J > l.QF due to finite order 
errors. 



D. Equal-time structure factor 

The momentum-resolved equal-time structure factor 
contains more information so that it is another interest- 
ing quantity 

5""(0) = ^Ee-^'3(/-O(^^"^p). (40) 

For a single quasi-particle it is directly connected to 
the full DSF by Eq. (jH]) because there is no mixing be- 
tween different quasi-particle spaces^. Our first focus 
is Sf^{Q). Within the CUT framework, this quantity 
can be computed from the effective observable cTj^g by 
Fourier transformation of the terms that exactly create 
one particle 

sriQ) = (0| <eff(-Q)<:Iff(Q) |0) (41) 



where n stands for the number of quasi-particles involved 
and a may take the values x or y. In Fig. [S] we compare 
the CUT results to the exact expression ([5a|) for the pa- 
rameters J = F, J = 1.5F and J = 1.9F. The agreement 
is very impressive though it worsens upon approaching 
the critical point. For J = F the DSF is essentially con- 
verged and the absolute errors are below 10~^''F"^. For 
J = 1.9F the error is below IQ-^F"! for Q < n/2. For 
Q > tt/2 the absolute error rises up to lO^^F^^. 

A closer analysis reveals that the largest absolute error 
occurs for all parameters at the wave vector Q = ir. The 
DSF diverges at this point for J ^ 2F. The relative 
error (not shown in the graphs) remains fairly constant 
over the whole Brillouin zone. Consequently our results 
differ from the exact ones found by Hamer et ali^ by only 
about 1% even close to the critical point at J = 1.9F. 

Finally, we consider S\^{Q). In Fig. |6]we compare the 
CUT results to the exact expression for the parameters 
J ^T, J = 1.5F and J = 1.9F. For J = F the DSF 
is essentially converged and the absolute errors are be- 
low 10~^''F~"'^. This changes for rising parameter J. For 
J = 1.9F the error is below IQ-^r-i for Q < 7r/2. For 
Q > tt/2 the absolute error remains below 10~'^F~^. In 
contrast to the S"^^ channel, the lowest absolute error is 
located in the channel at Q « tt. This is can be eas- 
ily understood because Sf^{Q ~ tt) constitutes a local 
minimum for all parameters J. Again, the relative error 
(not shown) remains essentially constant over the whole 
Brillouin zone. 



VI. DYNAMIC PROPERTIES 

In this section, we evaluate and present the dynamic re- 
sults for the transverse field Ising model. Here, 'dynamic' 
refers to frequency dependent quantities. We deal with 
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Figure 5. (Color Online) One-particle equal time structure 
factor Si^{Q) for the parameters J = T, J — 1.5T and J = 
1.9r. Comparison of the exact expression (jSap for the one- 
particle weight with the CUT results. 



Figure 7. (Color online) Energy dispersion for J — 2. OF (top) 
and J = l.QF (bottom). Comparison of the the CUT calcu- 
lations to the exact results. 
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Figure 6. (Color online) One-particle equal time structure 
factor S\^{Q) for the parameters J = T, J = 1.5F and J = 
l.QF. Comparison of the exact expression (|8b|l with the CUT 
results. 



the dispersion in Sect. IVI XI and with the DSF in general 
in Sect. IVI Bl and its different channels in Sects. IVI CI 
(5'^n. lVrD] (5'^^) and|VlE](S'««). 

The general DSF is an important quantity because it 
is directly measurable in scattering experiments. Fur- 
thermore, dynamic correlations strongly depend on the 
model under study and often exhibit features which re- 
veal the microscopic interactions in the Hamiltonian. De- 
spite the fact that the TFIM is integrable, the calculation 
of dynamic correlations remains a difficult and complex 
problem. 



A. Dispersion 



As before we were able to reach order 256 for the CUT 
calculation of the Hamiltonian. In particular, we obtain 
the hopping matrix elements up to a range of 256. For 
small parameters J, a low order calculation is sufficient to 
achieve a good agreement with the exact result. Closer to 
the critical point this changes distinctly, see Fig. [3 which 
makes higher order calculations necessary. For J ~ l.OF 
the absolute error of the result in order 6 is below O.OIF 
for q < 7r/2 and below O.IF for q > 7r/2. For the order 
32 result, it is below lO^^F for q < n/2 and below lO^^F 
for q > 7r/2. For J = 2r, the error of the order 32 result 
is below 10~3r for q < 7r/2 and below lO^^F for q > 7r/2. 
For the order 256 result, it is below lO^^F for q < Tr/2 
and below 10~^r for q > n/2. 

This behavior is expected because the excitations be- 
come more and more dispersive on increasing J. Con- 
sequently, hopping processes over larger and larger dis- 
tances become more important. To include these physical 
processes we need higher orders because the maximum 
range we can describe directly corresponds to the order 
of calculation (for lattice constant equal to unity). Di- 
rectly at the critical point the energy gap closes and the 
correlation length diverges concomitantly. 

The calculation of the dispersion is worst in the vicinity 
of the critical wave vector q = tt. We stress, however, 
that the value directly at g = tt, the energy gap of the 
TFIM, is calculated exactly up to numerical errors below 
j^g-iop rpj^-g jg accidental result because the energy 
gap happens to be a linear function of J so that it is 
captured correctly by any perturbative approach in linear 
order and beyond, compare Eq. ([?]). 
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B. Dynamic structure factor 

The DSF at T = is linked to the imaginary part of the 
retarded Green function by the fluctuation-dissipation 
theorem at zero temperature^ 



1 



ImG""(w,Q). 



(42) 



At T = 0, it is useful to write this Green function for 
a; > as a resolvent 



G""(c.,Q)==(g|a"(-g) 
1 



uj - {H{Q) - Eo) + lO 
where Eq is the ground state energy and 



a"(Q)|5) (43) 



-y 



(44) 



is the Fourier transformed spin operator cr". 

In the CUT framework, we replace all operators by 
the effective operators and the ground state by the zero- 
particle state, i.e., the vacuum of quasi-particles 



G""(L.,Q) = (0|5,"ff(-g) 
1 



{H,ff{Q)-Eo)+iO- 



-5Sf.(Q)|0). (45) 



We evaluate the resolvent in Eq. (|45j) by means of a Lanc- 
zos tridiagonalization yielding a continued fraction rep- 
resentation of the rcsolven t^^'^° 



G""(w,Q) 



(46) 




Figure 8. Qualitative sketch of the band-edge singularities in 
the DSF. 



to the continued fraction coefficients. 

For two massive hardcore particles without interaction 
and with finite range hopping in one dimension the band- 
edge singularities are known tobeQ; = /3 = l/2, see for 
instance Refs. [5l| and H^- We expect this behavior also 
to be true in the case of the TFIM because there is no 
interaction in the exact solution. In this case the relations 
P7a|) and (|47b)) yield 



1 



(49a) 
(49b) 



We confirm this behavior in the two-particle case of the 
S"^^ channel in Sect. lyTCl 



where the coefficients a„ and 6„ are the matrix elements 
of the tridiagonal matrix of the effective Hamiltonian. 
We refer the reader to Appendices [C] and |D] where we 
explicitly calculate S"ff{Q)\0) as well as the action of 
the effective Hamiltonian for the Lanczos tridiagonaliza- 
tion. The continued fraction is terminated by a standard 
square-root terminator. This is appropriate for square 
root singularities at the band-edges. 

Another piece of information that can be obtained from 
the sequences {a„} and are the exponents a and (3 of 
the band-edge singularities, see Fig. [51 They are directly 
connected to the asymptotics^ 



a-n = floo + Ooo h O I — 



(47a) 



br. = b^ + — + ^ j • ^47b) 

These relations allow us to obtain the band-edge singu- 
larities up to their signs by fitting 



fin) 



G+4 



(48) 



C. S^^ channel 

Because the observable stays local in the Jordan- 
Wigner representation of the TFIM the DSF in the S^^ 
channel can be obtained analytically, see Eq. ((7a)) . The 
DSF in the zz channel results from the two-particle con- 
tinuum. Even for large parameters J, no weight is shifted 
towards four or more particle spaces. All dynamics in- 
duced by this observable is captured in the two-particle 
sector. 

We emphasize that this fact holds as well in in the 
CUT treatment formulated in terms of the string opera- 
tor algebra. The corresponding local operator cr| = Oj^ 
is element of the string algebra so that its effective ob- 
servable after the CUT consists of a linear combination 
of string operators 

^i,cff = o-i+dO-i+d.o (50) 

d 

+ E _„ + h.c.) + 0++ ,„(0++ ,„ + h.c.] ) 
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where the maximum range n is hmited by the order of the 
calculation. We stress again that the operators Oj+d,o 
and O^^j^ „ do not create any excitations, while the op- 
erators O^^j^ „ create exactly two excitations. Thus, the 
vector S'"jj((5)|0) is only element of the zero- and of 
the two-particle Hilbert space. The same holds in the 
fermionic picture, where the operator is a local den- 
sity term which at most creates two fermionic excitations 
after the Bogoliubov diagonalization. 

Concomitantly, very high orders can be reached also 
in the transformation of the local observable. Because 
we transform a non translational-invariant operator, we 
have to consider the positions j + d and the starting site j 
so that the number of terms increases quadratically with 
the order. But we are still able to achieve an order of 
128. 

False color plots of the DSF obtained in this way are 
shown in Fig. IH] in order 128. The two-particle con- 
tinuum is depicted in dependence of the total momen- 
tum Q and the energy w. The overall intensity rises for 
larger parameters J . This stems from the decrease of the 
transverse magnetization which induces a shift of spectral 
weight from the zero-particle channel to the two-particle 
channel. Furthermore, we see that for small parameters 
J/(2r) most of the weight is concentrated in the region 
Q < tt/2 while the opposite behavior occurs for larger 
parameters x. Note the singularity inside the continuum 
on the right side of the Brillouin zone that separates two 
regions with low and high spectral weight, see the case 
J = 1.9r. Knowing the exact expression (|7a|) we can 
explain this singularity as van Hove singularity in the 
two-particle density-of-states. The two-particle energy 
uj[Q 12 -\- q) + LoiyQ /2 — q) displays a local maximum be- 
sides the global extrema as function of (/, if Q and J are 
large enough. 

We want to investigate more profoundly how the CUT 
calculation differs from the exact calculation by examin- 
ing the DSF for fixed parameters J and total momentum 
Q. Fig. [TU] shows S'^^(a;,Q) for J = 1.5r and J = l.QF 
and for the momenta Q ^ Q. Q — n/2 and Q ~ tt calcu- 
lated by the CUT in order 128 in comparison to the exact 
result. Note the excellent agreement for all parameters 
and momenta. The form of the DSF is very close to a 
half ellipse for low values of J because the continued frac- 
tion coefficients converge very quickly towards their final 
values Ooo and bca- For large J more spectral weight is 
concentrated at the lower band-edge which we attribute 
to a complex interplay between momentum and energy 
conservation. For the parameters J = 1.9 and Q = 7r/2, 
one clearly sees the singularity inside the continuum of 
the DSF which is the above mentioned van Hove singu- 
larity from a local maximum. 

A detailed analysis shows that the error is lower in the 
middle of the continuum than at the band-edge singulari- 
ties. This is expected due to the strong change of the DSF 
at the edges. On average, the error is below 10~^r~^ 
even for large parameters J. At the band-edges the error 
can rise up to 10~^r^"'^. We presume that the errors are 
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Figure 9. (Color online) Tlie DSF S'^^((.j,Q) for the parame- 
ters J = r (top), J = 1.5r (center) and J = 1.9r (bottom). 
The maximum range for the Lanczos algorithm is dmax = 1000 
sites, the continued fraction was evaluated to a depth of 50 
and then terminated by the square root terminator. The color 
indicates the spectral density, see legend to the right. The up- 
per and lower edge of the two-particle continuum are indicated 
by white lines. 



mainly due to inaccuracies in the Lanczos tridiagonaliza- 
tion and due to the limited maximum range in the trans- 
formation of the observable by the CUT. Nonetheless the 
errors are still very small and justify our approach. 

Next, we investigate how the continued fraction coef- 
ficients approach their limiting values. Thereby, we es- 
timate the exponents of the band-edge singularities ac- 
cording to Eqs. (|47ap and (j47bp . The continued fraction 
coefficients for the case J = l.SF and total momenta 
Q = and Q = 7r/2 arc shown in Fig. [11] The coeffi- 
cients for the case Q — approach their limit exponen- 
tially. Therefore we know by Eqs. (|47ap and (j47bl) that 
both exponents take the value 1/2. 
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Figure 10. (Color online) DSF S^^{iL},Q) for the parameters 
J = 1.5r (left) and J = 1.9r (right) for three total momenta 
Q. The maximum range for the Lanczos algorithm is dmax = 
4000 sites, the continued fraction was evaluated to a depth of 
100 and then terminated by the square root terminator. 
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Figure 11. (Color online) Absolute difference between the 
continued fraction coefficients and their final values for the 
case J — 1.5T and total momenta Q = and Q = n/2. 



For the case Q = 7r/2, the coefficients do not converge 
so rapidly. We fit them for this case versus 1 /n^ to check 
if they display terms in 0{l/n^). Both coefficients oscil- 
late around the final value which can not be described 
by Eqs. (|47ap and (j47b[) . This again indicates that the 
exponents at the band-edges are 1/2. We also checked 
other momenta and they support the assumption that all 
exponents are 1/2 for the two-particle case as it has to 
be according to the fermionic analytical results. Thus, 
these findings corroborate the validity of our approach 
and analysis. 



D. 5"==" channel 

In Ref. [s^ Hamer et al. derived an analytic expres- 
sion for the one-particle contribution for the longitudinal 
structure factors. To our knowledge no data is available 
in the literature for higher quasi-particle contributions 



Q 


a 






7r/2 

TV 


2.5 ± 0.3 
3.0 ± 0.2 
3.0 ± 0.2 


1.0 ± 0.2 

2.7 ± 0.2 

2.8 ± 0.1 



Table II. Exponents for the band-edge singularities of 
S3^{iJ,Q) for J = l.Sr. The errors are determined from the 
fits using the Levenberg-Marquardt algorithn^'^ 



away from the scaling rcgior^. Here our approach is 
able to provide complementary quantitative knowledge. 

Similar to the two-particle case S^^{uj,Q), the three- 
particle case S^'^{uj, Q) consists of a continuum of states. 
We are limited to a maximum order 38 due to the com- 
plicated structure of the local observable which is not 
part of the closed string algebra. Overview plots for the 
DSF obtained by the CUT are found in Fig. HH In these 
plots the three-particle continuum is depicted in depen- 
dence of total momentum Q and the energy lo. 

The total weight rises on increasing J because spec- 
tral weight is transferred from the one-particle sector to 
the higher quasi-particle channels. The same qualitative 
behavior is observed for dimerized spin chains and spin 
ladders, and related systemsS^i^. In addition, we no- 
tice that most of the spectral weight is concentrated at 
momenta Q < 7r/2 for small parameters J. The weight 
slowly shifts for growing J similar to the case. For 
J = 1.9r, most of the spectral weight is concentrated 
rather strongly at the lower band-edge of the continuum. 
The same tendency was found in the S*^^ case as well. 
Still the shape of the DSF in the 5"^^ case differs strongly 
from a semi-ellipse in contrast to the DSF. 

A more quantitative investigation is shown in Fig. 1131 
where S^^ is plotted for J = l.QF and momenta Q ~ 
0, Q = 7r/2 and Q ^ it. It is confirmed that most of 
the spectral weight is concentrated at the lower band- 
edge for large values of J. For Q = tt, a strong wiggling 
occurs which is to be attributed to the errors due to the 
calculation in finite order. 

Next, we investigate the band-edge singularities in the 
three-particle case S^^{uj, Q). As in the two-particle case, 
we use the relations (|47ap and (|47ap by fitting a 
power law to the continued fraction coefficients. An ex- 
ample is shown in Fig. [TJ] 

In contrast to the two-particle case, no exponential ap- 
proach towards the limit values occurs. For all momenta, 
both a„ and 6„ show a behavior proportional to l/n^. 
We stress that the 0{l/n^) terms are significant up to 
l/n"^ « 0.0002 ^ n « 70. The values for the band-edge 
singularities obtained from the fits are shown in Tab. HIl 

For Q = TT and Q = 7r/2, both exponents are close to 
3 while for Q = the exponents differ and we deduce 
that a = 2.5 and /3 = 1 holds. We stress, however, that 
the obtained exponents may still be affected by rather 
large numerical errors. In Ref. Issl a general expression 
for the multi-particle band-edge singularities is derived 
for a simple one-dimensional model of hardcore bosons 
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Figure 12. (Color online) The DSF 5'f"^(a;,Q) for the param- 
eters J — r (top), J — 1.5T (center) and J — 1.9T (bottom). 
The maximum range for the Lanczos algorithm is dmax = 100 
sites, the continued fraction was evaluated to a depth of 50 
and then terminated by the square root terminator. The color 
indicates the spectral density, see legend to the right. The dis- 
persion is indicated by the white solid line. The upper and 
lower edge of the three-particle continuum are indicated by 
white dashed lines. All results are computed in order 38. 
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Figure 13. (Color online) DSF S§^{uj,Q) for the parameter 
J = 1.9r for three chosen total momenta Q. The maximum 
range for the Lanczos algorithm is dmax ~ 300 sites, the con- 
tinued fraction was evaluated to a depth of 100 and then ter- 
minated by the square root terminator. 
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Figure 14. (Color online) Continued fraction coefficients for 
the case J = l.SP and total momentum Q — 0. The upper 
panel shows the coefficients On and the lower panel shows the 
coefficients 6„. The limit values are indicated by horizontal 
lines. The red/green lines indicated linear fits in 1/n^; note 
the scale of the a;-axis. 



with nearest-neighbor hopping. The result reads 

Sn oc oj 2 for n > 1, (51) 

, but close to the extrema of the dispersion. For the 
three-particle case, this yields Sn oc which agrees well 
with our results for Q = tt and Q = n/2, but differs 
for (5 = 0. The discrepancy in the latter case may be 
due to the more complex structure of the dispersion in 
the effective Hamiltonian for the TFIM which includes 
longer range hopping processes. 



E. 5"^^ channel 

As in the S'^^ channel, the S^^ channel consists of 
1,3,5,... particle contributions. Overview plots for the 
three-particle DSF obtained by the CUT in order 38 are 
found in Fig. [151 In these plots, the three-particle con- 
tinuum is plotted in dependence of the total momentum 
Q and the energy uj. For small J, the S^^ channel looks 
similar to the channel. The only difference is the ab- 
solute weight because the three-particle continuum in the 
S"^"^ channel gains weight sooner, i.e., for smaller J/{2T), 
than in the S^^ channel. 
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Figure 15. (Color online) DSF S^^{u],Q) for the parameters 
J = r (top), J = l.Sr (center) and J = l.QF (bottom). The 
maximum range for the Lanczos algorithm is dmax = 100 sites, 
the continued fraction was evaluated to a depth of 50 and then 
terminated by the square root terminator. The color indicates 
the spectral density, see legend to the right. The dispersion is 
indicated by the white solid line. The upper and lower edge 
of the three-particle continuum are indicated by white dashed 
lines. All results are computed in order 38. 



For higher values of J, there are already qualitative 
differences between the S^^ and the S'^'^ channel. Most 
of the spectral weight is still concentrated at the lower 
band-edge of the continuum. No spectral weight is gained 
in the region of the critical wave vector Q = tt which 
constitutes a major difference to the S^^ channel, see 
Fig.HH 

Scans of 5*1^ at fixed Q are shown in Fig. [1^ for J = 
1.9r and momenta Q = 0, Q = tt/2 and Q = tt. For 
Q < 7''/2, most of the spectral weight is concentrated at 
the lower band-edge. This changes distinctly for Q > 
7r/2, especially for Q k tt. Here spectral weight is spread 
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Figure 16. (Color online) DSF S^^{uj,Q) for the parameter 
J = l.QF for three chosen total momenta Q. The maximum 
range for the Lanczos algorithm is dmax ~ 300 sites, the con- 
tinued fraction was evaluated to a depth of 100 and then ter- 
minated by the square root terminator. 



Q 


a 


P 




7r/2 

TT 


2.9 ± 0.1 
2.9 ± 0.2 
2.9 ± 0.1 


1.0 ± 0.1 

2.7 ± 0.2 

2.8 ± 0.2 



Table III. Exponents for the Band Edge Singularities for 
S^^(u],Q) for J = 1.5r. The errors are determined from the 
fits using the Levenberg-Marquardt algorith m^^'^^ . 

rather equally over frequency space. We also observe 
some wiggling which is due to finite order errors. 

The values for the band-edge singularities obtained by 
fits as described before in the channels S^" and Sf^ are 
given in Tab. lIIll They mostly equal those obtained for in 
the channel within numerical errors. Only the case 
(5 = differs. We deduce that a = 3 and 13 = \ holds for 
S^^ generally. 

VII. CONCLUSION 

Summarizing, we showed that the one-dimensional 
transverse field Ising model (ID TFIM) in the high 
field limit can be expressed in terms of string operators 
which form an algebra which is closed under commu- 
tation, which agrees with the previous finding by Jha 
and Valatir^. This property allowed us to solve the 
ID TFIM in the high field limit to very high accuracy 
by continuous unitary transformations without resorting 
to the Jordan- Wigncr transformation to non-interacting 
fcrmions. Note that the solution provided formally also 
covers the low field limit due to the duality of the model. 

The only remaining restriction in the presented solu- 
tion is the truncation in a given order n in the ratio 
X = J/(2r) of the Ising coupling J to the field strength 
r. But due to the closure of the string algebra the num- 
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ber of terms to be tracked grows only linearly in the order 
n so that very high orders up to n = 256 can be achieved. 
Thus, accurate results for all practical purposes could be 
obtained. The order corresponds directly to the range of 
physical processes which are included if the lattice spac- 
ing is set to unity. We employed the recently developed 
deepCUT approach which is pcrturbatively correct in the 
targeted order and provides a robust extrapolation be- 
yond this order—. 

High orders are accessible for the Hamiltonian and all 
observables which belong to the string algebra. They 
cannot be obtained for observables which do not belong 
to the string algebra such as the longitudinal spin com- 
ponents. For this reason, the longitudinal components 
could be unitarily transformed only up to order 38. 

We gauged the results in computing various static 
properties such as the ground state energy, the trans- 
verse magnetization, and the one-particle contribution 
to the equal time structure factors Si"{Q) for a = x 
and a ~ y. The first two quantities can be compared 
directly to the analytically accessible results via Jordan- 
Wigncr transformation. The onc-particlc contribution to 
the equal time structure factors has been conjectured by 
Hamer and co-workers by series expansions and strongly 
underlined by the mapping to a two-dimensional classical 
Ising model for which the exact results are known^^. 

Similarly, we computed dynamic quantities such as 
the one-particle dispersion and the momentum- and 
frequency-resolved diagonal dynamic structure factors. 
The dispersion and the transverse structure factor can 
again be gauged against the analytical result obtained 
in terms of non-interacting fermions. The longitudinal 
structure factors are much more difficult to address be- 
cause their excitation operators are highly non-local in 
terms of the non-interacting fermions. While the one- 
particle contribution can be derived from the static one- 
particle structure factor and the exactly known disper- 
sion there arc no results for the next important three- 
particle contributions for general coupling x' < 1. Only 
in the scaling region around the quantum phase tran- 
sition results for the three-particle contributions in fre- 
quency domain exisli^. Our results are reliable further 
away from the scaling region so that they are complemen- 
tary to the existing information. The equation of motion 
approach pursued by Perk and Au-Yangii provides in- 
formation on the correlations in the time and real space 
domain. But so far no analysis with frequency and mo- 
mentum resolution has been performed. 

The presented theoretical thrce-particlc data for the 
static and the dynamic structure factor provides predic- 
tions where in momentum and frequency space one can 
expect significant three-particle signal. This information 
may guide future experimental searches for many-particle 
contributions. Concretely, our results show that the S^^ 
channel is considerably better suited for such searches 
than the S^^ channel. In the S^^ channel the single 
particle contributions dominates over the multi-particle 
contributions except very close to the quantum phase 



transition. 

Moreover, we found that the spectral weight in the 
three-particle dynamic structure factors is concentrated 
close to the lower band-edge if the parameters are such 
that the system is not too far away from criticality. Fur- 
ther away from criticality the main response is rather 
featureless and hardly displays a dependence on the total 
momentum Q. Then the spectral weight is still concen- 
trated close to the lower bandedge around (5 = 0, while 
it is spread out in the middle of the band around Q = tt. 

Our approach can be pursued further for all one- 
dimensional models of which the Hamilton operators can 
be expressed within the string algebra. Further investi- 
gations for other response functions are possible as well. 
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Appendix A: Closure of the string algebra 

Here we show that the string algebra is closed under 
commutation. This means that the commutator of two 
string operators can again be written as a linear com- 
bination of string operators. First, we show that on a 
chain two string operators commute if neither of their 
start-/cnd-opcrators arc on the same site. Without loss 
of generality this means 



= 
= 
= 



j,n' ^l.rn 
r)<t"^ QXi 



I > j, I + m < j + n, (Ala) 
I > j, I + m > j + n, (Alb) 
l> j + n, (Ale) 



for (j), e, X, C S {+, — }. The last commutator (jAlc|) is zero 
because operators acting on completely different sites al- 
ways commute in a bosonic algebra. A simple calculation 
yields for the first commutator (jAlap 



Q-P^ QXi 



'I ^l+mi"l 



(A2a) 
(A2b) 



(A2c) 

I, (A2d) 
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using [(7^, a^] = y^a^^ a'^a^ = xcr-^ and a^a^ = — xcr-^, if 
all operators act on the same site. Analogous calculations 
yield that (jAlb|) holds as well. 

The remaining contributions consist of commutators 
where either the start- and/or end-operators are on the 
same site. The start- and/or end-operators on the same 
site of two string operators need to be different because 
otherwise the identities a'^a'^ = a~a~ = imply a 
vanishing result. Explicit calculations yield for the non- 
vanishing commutators 



3,n ' j,m 



- fO^' 

S J +771, n—m' 



with m < n and 



j,7i+m' 



n,0 



~e20' 



Q0£ Q-4>- 



(A3) 

(A4a) 
(A4b) 
(A4c) 



Explicit calculations for the translationally invariant 
string operators yield the following set of commutator 
relations for the case n,m € N^, n < m 





T — 

rn 


_|_ h — 

n+m ' n+'/?i m—n 


^771^71' (A5a) 




T+-' 


_ _ 7^ + + 


(A5b) 




T-+' 

rrt- 


_ ^ + + _ + 

n+m rn — n ' 


(A5c) 




^ n 


T^ + + 1 T^+ + 

n+m ^ 711 — 711 


(A5d) 






T^ + + 1 T^+ + 

n+m "T" 771 — ni 


(A5e) 




To ] 




(A5f) 




771 


-0, 


(A5g) 



and for the case 



\ rn ' 777 



2m 



Tr 



2m 



-To, 



L 777 5 777 
[71 + + rp- + 

I m ' 777 

[7^-1-+ 7^++ 

L m ' ^ 777 



71+-!- 

^2m I 

0, 



(A6a) 
(A6b) 
(A6c) 
(A6d) 



which are all linear combinations of string operators. 
Contributions with n > m are also included by exchange 
of the arguments in the commutators. This concludes 
the derivation of the closure of the string algebra. 



of maximum range n, Tq, T^^,T^^, . . . we can cre- 
ate a new complete set of operators of range n -I- 1 by 
commutation with a string pair-creation-operator, 



1 _ rj.+ - 



= r 

71 I J 

_ 7^+ + 

^ -^71 + ] 

= -T" 



(Bla) 
(Bib) 

+ T-I,. (Blc) 



Thereby, we created the string operators of range 71 + 1. 
Because the Hamiltonian in Eq. (j3ip already comprises a 
complete set of range one we can deduce that all ranges 
n G N+ will be created during the flow. Hence, we can 
conclude for the generator of the TFIM 



5:4+ (r++-T--) 



(B2) 



For step two we consider the relations in Eq. (|A5|) and Eq. 
(|A6[) . We start with the contributions to the operator Tq. 
Such contributions are created only in the case m = n. 
For a given range n there are two contributions from the 
commutators 



\t++ t — 

[ 77 1 71 

r ++n 
\_-^ n 1-^71 J 



= To 
= -To 



(B3a) 
(B3b) 



both with prefactor one. Note that T^^ and have 
the same prefactor up to a sign due to hermiticity/anti- 
hermiticity. Finally, these considerations yield 



dit 



71=1 



+ + \ 



(B4) 



for the flow equation for the prefactor of Tq. 

Next; we consider the operator T^~ and T^^, respec- 
tively. They are created by two different kinds of com- 
mutators. For k + I = m 



and for \k ~ l\ = m 



= T+- 

7n 



[Tft^,Tr] = -T+- 



(B5a) 
(B5b) 

(B6a) 
(B6b) 



with prefactor one. Note that the operator T^"*" have the 
same prefactor as T^~ . These calculations yield 



Appendix B: Proof of infinite order flow equation 

To prove the expression ([55)) we proceed in two steps. 
First, we show that all kinds of string operators of ar- 
bitrary range will be created during the flow. Next we 
show which contributions occur in the DES. 

Our starting point for step one is the Hamiltonian of 
the TFIM in string operators in Eq. ([51]) . By induction 
we show that once we have a complete set of operators 



k-\-l—7n 
k,l 



k 



\k — l\ —m 
k,l 



(B7) 



Last we consider the operator and T",^^, respec- 

tively. They are created by three different kinds of com- 
mutators. For k + I = m 



[T^\T^ 



77- 



(B8a) 
(B8b) 
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for \k — l\ = m 

[r++, r+-] = sgn(fc - l)T++ + ..., (B9a) 
[T++, T-+] = sgn(fc - l)T++ + . . . , (B9b) 

where the sign function stems from the different signs 
in the cases [T++,T+-] and [T++,T+-] in Eq. IMi- 
Finally, the third case is given by 

[T++,To] =-4T++. (BIO) 

Now we can write down the complete flow equation for 
the prefactor t++ 

k-\-l—m 

k,l 

\k—l\—m 

+ 2 sgn{k-l)t++t+-, (Blla) 

k,l 

which concludes our derivation for the flow equation for 
infinite order. 



Appendix C: Calculation of S"g(Q)\0) 

We start from Eq. (j45|). To apply the Lanczos algo- 
rithm we need to calculate 

%(Q)|0). (CI) 

We split the vector into its components of different par- 
ticle number. For the two-particle structure factor the 



state 

S^,\liQ)\0)^^ E ^'""^ 

r,do,di,j 

|r + do,r + rfo + rfi), (C2) 

with di > must be considered. The sum over j ad- 
dresses all operators that create an excitation at r + 
and another at r + da + di which are different in their 
content of factors erf at various sites. The index j is used 
to distinguish them. In contrast, in a strict multi-particle 
representation there would be only one operator. Shift- 
ing the exponent by dQ + di/2, the center of mass, results 
in the expression 

do,di,j 



:=|Q.rfl> 

(C3a) 

= s,^«;f(Q)|g,di), (C3c) 

where we have introduced \Q,di) which is the Fourier 
transformation of a two-particle state with distance di . 
For the three-particle structure factor the state 

5eff|^(Q)|0)=-^ Y ^'^'-^^'^ 
7\dQ.,di ,d2,j 

■ \r + do, r + (io + di,r -\- do di -\- ^2) , 

(C4) 

with di > and d2 > must be considered. Shifting 
the exponent by do + 2di/3 + d2/3, the center of mass, 
results in the expression 



5'cff|o(Q)|0) = J2 e-*«W-'+2'^^/3+'^^/3)^do.di.d2 



do,di,d2,j 



^ Y e^Qir+d„+2d^/3+d2/3} \r + do,r + do + di,r + do + di + ds), (C5a) 



■. = \Q.di,d2) 

J2 e-*«(*+2'^i/3+'^^/3)^do.jiA |g^^^^^^^^ (C5b) 
do,di.d2 ,j ^ 



Y st"gf'''^Q)\Q,d,,d2), (C5c) 

do,di,d2,j 
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where we introduced |Q,(ii,d2), which is the Fourier 
transformation of a three-particle state with distance di 
between the first two particles and distance ^2 between 
the second two particles. 



Appendix D: Action of the effective Hamiltonian 

To apply the Lanczos algorithm we need to know the 
action of the effective Hamiltonian on the two- and three- 
particle states, calculated in App.[C] We stress that after 
the CUT there are no terms that violate particle-number 
conservation. We analyze the action of the operators 
Tq^T^ ,T~'^ separately. Starting with the simple oper- 
ator To yields 



I 

io(oo)To \Q, di) = to(oo)ro^ ^ g.Q(r+do+di/2) |^ + ^q, r + do + rfi) , (Dla) 

V ^ 

= to(oo)(-iV + 4)^ e*^(''+*'+'^^/2) \r + d„,r + do + rfi) , (Dlb) 
^ iEo + Uo{(X)))\Q,di) . (Die) 

I 

Next, we analyze the action of the operator 

I 

^t+-(oo)r+- |Q,di) ^ ^(-l)"-ii+-(oo)e^'^"/' |Q,di +n) 

n n 

+ E (-l)""'C"(oo)e*'5>V2 |Q,di -n) (D2a) 

n<di 

I 

and of the operator T~~^ 

I 

E t+-(oo)r-+ IQ, di) = E(-l)"-4+-(oo)e-'0"/2 \Q, d, + n) 

n n 

+ J2 (-l)""'i;l"(oo)e-'«"/2 \Q,di - n) (D3a) 
+ 5](-l)%t"(oo)e-*'3"/2|Q,n-di). 

n>di 

I 

Note the different signs of the second and third term due Next, we analyze the action of the operator 
to the properties of the string operator. 

Similarly to the two-particle state we examine the ac- 
tion of the effective Hamiltonian on the three-particle 
state. The simple operator To yields 

io(oo)To IQ, di, ds) = to((X3)(-iV + 6) |Q, di, da) , 

(D4a) 

= (i;o + 6to(c»))|Q,di,d2). (D4b) 
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^ <.r (c»)r+- IQ, d,,d2) = 5^(-l)"-4+-(c»)e'Q"/3 \Q, di + n, ds) 

+ H (-l)""'4"(oo)e*«"/3|g^^^_„^^^^^^ 

n<(/i 

+ ^(-irt+"(oo)e''^"/'IQ,^-rfi-'^2 + di) (D5) 

n>di 
n<d2 

+ E (-ir4"(oo)e*'^"/3|Q,di+(i2-n,n-(i2)) 

di+d2>n>d2 
n>di+d2 

I 

and of the operator Tj7+ 

I 

tt-{^)T-+ \Q, di, d2) = 5^(-l)"-i<+-(oo)e-''?"/3 \Q, d^,d2+ n) 

n n 

+ E (-l)""'4"(oo)e-^«"/3 IQ, di +n,d2~ n) 

n<d2 

+ (oo)e-'Q"/3 IQ, di + d2, n - ^2) (D6) 

n>d2 

+ E (-l)""'4"(oo)e-'«"/3 IQ, di - n, ^2) 

+ Y {-irt+-{°o)e-"^''^^\Q,n-di,di+d2-n) 

di+d2>n>di 

+ E (-ir"'4"Me-^'5"/3|Q,d2,n-di-d2>. 

n>di+d2 

I 
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